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We combine computer simulations and analytical theory to investigate the glassy dynamics in 
dense assemblies of athermal particles evolving under the sole influence of self-propulsion. Our sim¬ 
ulations reveal that when the persistence time of the self-propulsion is increased, the local structure 
becomes more pronounced whereas the long-time dynamics first accelerates and then slows down. 

We explain these surprising findings by constructing a nonequilibrium microscopic theory which 
yields nontrivial predictions for the glassy dynamics. These predictions are in qualitative agreement 
with the simulations and reveal the importance of steady state correlations of the local velocities to 
the nonequilibrium dynamics of dense self-propelled particles. 
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I. INTRODUCTION 

The application of statistical mechanical methods to 
the dynamics of individual motile objects started shortly 
[I] after Einstein’s work on Brownian motion [2]. Re¬ 
cently, the collective behavior of systems consisting of in¬ 
teracting self-propelled particles attracted interest 13 ®- 
An important motivation for studying ‘active’ systems is 
to understand spectacular dynamics observed in assem¬ 
blies of living systems, such as coherent motion [30. A 
more fundamental motivation stems from the nonequi¬ 
librium nature of active systems that are driven by inter¬ 
nal, non-thermal self-propulsion forces, which represents 
a difficult challenge for statistical physics. The behavior 
of such systems may defy our equilibrium-based phys¬ 
ical intuition, as demonstrated by large-scale collective 
motion in persistent hard disk systems 7], and the emer¬ 
gence of dynamic clustering Eli, phase separation 0- 
m, and nonequilibrium equation of state HMD in re¬ 
pulsive self-propelled particles. 

Active particles may also exhibit behavior similar to 
that found in equilibrium systems, such as crystallization 
m ■ This behavior suggests that, like thermal systems 
EH, dense active systems could possibly be supercooled 
and exhibit glassy dynamics; as recently argued theoret¬ 
ically on the basis of a simple driven glassy model P3- 
This study was followed by numerical investigations of ac¬ 
tive particles with hard-core mang or continuous [2DU2I] 
interactions. Active glassy dynamics was observed, but 
when compared with thermal systems, the onset of glassy 
behavior was always pushed towards higher densities and 
lower temperatures. When self-propulsion is progres¬ 
sively added to an otherwise thermal system, it was noted 
that the local structure becomes less pronounced [l8j|2Tj. 
Such change in local structure suggests both a simple 
physical explanation for the shift of the onset of glassy 
behavior, and that a straightforward extension of mode¬ 
coupling theory (22] can describe this shift 17, 23 . 

Here we present a simulational and theoretical study of 


the structure and glassy dynamics of a more complex sys¬ 
tem in which self-propulsion is the only source of motion. 
We study self-propelled particles interacting with a con¬ 
tinuous potential (differently from [l2J[]]I]) and without 
thermal Brownian motion (differently from jTB jfl8l 120] ). 
Thus, our model is ‘athermal ’ EHIM], and the degree 
of nonequilibrium is quantified by the persistence time 
t p of the self-propulsion. As t p increases at a constant 
effective temperature, we find that the structure of the 
system becomes more pronounced, whereas the dynam¬ 
ics initially speeds up and then slows down, showing that 
in our model the shift of the glassy dynamics with self¬ 
propulsion is not simply the direct consequence of the 
changing microstructure EQE3. 

To elucidate these findings we develop a microscopic 
theory that accounts for the nonequilibrium nature of 
athermal self-propelled systems. We are aware of no 
other approach where many-body interactions are taken 
into account at the microscopic level. The theory is con¬ 
structed from the steady state structure factor and the 
steady state correlations of the velocities. The latter cor¬ 
relations are non-trivial only for self-propelled systems 
with a finite persistence time, and are central to explain 
the opposite, and seemingly contradictory, evolution of 
the structure and dynamics in dense active materials. 

II. MODEL ACTIVE SYSTEM 

We study a system of interacting self-propelled parti¬ 
cles in a viscous medium. We model self-propulsion as an 
internal driving force evolving according to the Ornstein- 
Uhlenbeck J25J process: 

U = ^[Fi+f-], (1) 

fi = -Tp'ti +T]i. (2) 

In Eq. ([I]), r.i is the position of particle i , £o the fric¬ 
tion coefficient of an isolated particle, F* is the force 
acting on particle i originating from interactions, F; = 


2 


— an d i s the self-propulsion acting on 

particle i. In Eq. ([2]), r p is the persistence time of the self¬ 
propulsion and r/i is an internal Gaussian noise with zero 
mean and variance ( J 7i(t) I 7j(t , )) no i se = 2 DfISijS(t — t'), 
where (...} noise denotes averaging over the noise distri¬ 
bution, Df is the noise strength and I is the unit ten¬ 
sor. Without interactions, Eqs. (l]|2) produce a persis¬ 
tent random walk with a self-diffusion coefficient Z?o = 
H/Tp/^g, which defines the single-particle effective tem¬ 
perature: T e s = D 0 £ 0 = DfTp/t ;o [33] . It is convenient 
to choose as three independent control parameters the 
number density p (which is kept fixed), the effective tem¬ 
perature Teff, and the persistence time t p . The persis¬ 
tence time quantifies the degree of nonequilibrium; when 
t p 0 our system becomes equivalent to a Brownian 
system at temperature T = T e g. Including additional 
thermal noise would add a fourth control parameter to 
the model. Finally, we should mention that recently an 
approximate mapping has been derived m between our 
model and the standard active Brownian particles model 

ESI. 


III. COMPUTER SIMULATION STUDY 

To compare the glassy behavior of the self-propelled 
system with that of a well-studied thermal system, we 
simulated the Kob-Andersen (KA) binary mixture | 29 | . 
All quantities presented pertain to the large particles, 
which comprise 80% of the mixture. The results are pre¬ 
sented in reduced units [3UI at the well-studied density 
p = 1 . 2 . In panels (a-c) of Fig. [l]we show the dependence 
of the structure and dynamics when moving away from 
equilibrium by increasing the persistence time of the self¬ 
propulsion at constant T e ff = 0.5 [31]. At T = 0.5 the 
thermal KA system exhibits glassy dynamics. In panels 
(a,c), we see that the pair correlation g(r) of the active 
fluid becomes more structured at all length scales with 
increasing t p . In equilibrium, such behavior is usually 
accompanied by slower dynamics. Surprisingly, panels 
(b,c) show that the nonequilibrium dynamics exhibits a 
non-monotonic dependence on t p , which allows one to de¬ 
fine an ‘optimal’ value for t p . Describing the contrasting 
dependencies of structure and dynamics is a theoretical 
challenge since most microscopic glass theories predict 
the dynamics on the basis of the pair structure mm- 
A reentrant behavior of the dynamics in driven glassy dy¬ 
namics is typically not observed P2H2H- In panels (e,f) 
we show that, if the persistence time is chosen at its op¬ 
timal value, the dependence of the dynamics of the self- 
propelled system on T e g is weaker than that of the ther¬ 
mal Brownian system on T. In particular, we can study 
the self-propelled system at T e g = 0.4, whereas it is chal¬ 
lenging to simulate the thermal system below T ss 0.43. 
The opposite is true for longer persistence times, where 
the dependence of the dynamics of the self-propelled sys¬ 
tem on T e ff becomes significantly more pronounced than 
that of the thermal Brownian system. 




FIG. 1: Evolution of the structure and dynamics with the 
persistence time t v at T e g = 0.5 (a-c) and comparison of dy¬ 
namics in active systems with r p = 2 x 10 -3 and Brownian 
systems at T = T e g = 0.5 (d-f). In (a,b) solid lines correspond 
to t p = 0 (equivalent to a thermal system at T = 0.5), and 
the dashed, and dot-dashed lines correspond to t p = 2 x 10~ 3 , 
and 2 x 10 2 , respectively. In (a) we show the pair distribu¬ 
tion function g(r)- the curves are shifted vertically for clar¬ 
ity. In (b) we show the self-intermediate scattering function 
F s (q-,t) = (e* q ' < ' r A t f _r j(°))} for q = 7.25. In (c) we show the 
relaxation time, T a (circles) and the peak value of g(r ) (tri¬ 
angles). The relaxation time is defined as F s (q-,T a ) = e -1 . 
In (d,e) solid lines correspond to active systems at T e g = 1, 
0.6, 0.45 and 0.4, and dashed lines correspond to equilibrium 
systems at T = 1, 0.6 and 0.45 (from left to right). In (d) we 
show mean-square displacement (<5r 2 (f)) = ((r j(t) — rj(0)) 2 ) 
and in (e) we show F s (q ; t). In (f) circles represent r a , squares 
represent the inverse self-diffusion coefficient, 0.1 /D of active 
(closed symbols) and thermal (open symbols) systems. 


IV. THEORY 


We now outline a microscopic theory for the time de¬ 
pendence of the collective intermediate scattering func¬ 
tion, F(q ; t), of our model active system (331 . 

F(q;t) = N- 1 (n( q )e nt n(-q)) . (3) 

In Eq. (§, N is the number of particles, n(q) = 
yu e _lq rj is the Fourier transform of the microscopic 
density, and Q is the V-particle evolution operator that 
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can be derived from the equations of motion 00: 

« = v Ht + F <)+E | - (Vf< + £/|-) ■ 

i i x 7 

(4) 

Finally (...) in 0 denotes an average over a steady 
state distribution of positions and self-propulsions; the 
steady state distribution stands to the right of the quan¬ 
tity being averaged, and all operators act on it too. In 
our approach, we first integrate out the self-propulsions 
and then we use the projection operator method, and a 
mode-coupling-like approximation to derive an approxi¬ 
mate equation of motion for F(q; t). 

To begin, we briefly discuss the case of non-interacting 
particles. In this case one could start from the Laplace 
transform of Eq. © with the evolution operator similar 
to 0 but without interactions. In the simplest approxi¬ 
mation, after integration over self-propulsions one gets 

F(q;z)=N- 1 (n(q)(z-Q,g ee (z)) 1 n(- q)^, (5) 


where fl| f ee (z) = £ 0 2 V* (z + r p *) 1 D f r p • V* and 
(...} r denotes the steady state average over particles po¬ 
sitions. According to Hf® e (z), particle motion is ballistic 
at short times and diffusive at long times, with the long¬ 
time self-diffusion coefficient D 0 discussed above. 

For interacting particles, the integration over self¬ 
propulsions is more complicated due to non-trivial cor¬ 
relations between positions and self-propulsions (already 
present for a single self-propelled particle in an external 
potential [26]). As we show in Appendix A, the final 
result is a formula analogous to Eq. ([5]) but with the fol¬ 
lowing evolution operator, 

i,3 

(6) 

Here, (•■•) lss is an average over the conditional 
steady state distribution of self-propulsions, 
^’|f(ri,fi,...,r J v,fjv)/^(r 1 ,...,riv), where the su¬ 
perscript ‘ss’ stands for ‘steady state’. Furthermore, 
F° ff = Vj lnPjy(ri,..., rjv) is the effective force acting 
on particle j in the steady state. The most important 
physical assumption used in the derivation of Eq. © 
is the absence of systematic currents in the steady 
state, see Appendix A. We expect the persistence time 
to be renormalized by the interactions; the presence 
of the bare persistence time in Eq. © represents an 
approximation. 

We now use the projection operator method and arrive 
at the following memory function equation, 

d?F{q; t) + T~ 1 d t F(q; t) = (7) 

~ ^ls(g) g F ^ q; J dt ' M1 ~ ^ dt ' F ( q ’^)■ 


Here, S{q) = (n(q)n(—q)) is the steady state structure 
factor, wy (q) = (N£$) 1 <^|q • J2i ( f * + F ») e“ iq ' ri | 2 ) 

quantifies correlations of the velocities of individual par¬ 
ticles [32], and M 1TI (q; t) is the irreducible memory func¬ 
tion. The presence of the second time derivative in 
Eq. 0 reflects the ballistic nature of the short-time mo¬ 
tion. A comparison of Eq. 0 with the analogous equa¬ 
tion for an underdamped thermal system suggests to in¬ 
terpret t p u)\\ (q)/S(q) as a short-time collective diffusion 
coefficient. Since this coefficient involves W||(<?), even in 
the absence of the memory function we need two static 
correlation functions to predict the dynamics, S(q) and 
ai\\(q). Whereas the emergence of velocity correlations 
can be generically expected far from equilibrium, the spe¬ 
cific role they play for self-propelled systems is non-trivial 
and was not identified before. 

The main approximation of our theory is a factoriza¬ 
tion approximation for the memory function to close the 
dynamical equations. This is analogous to the mode¬ 
coupling approximation [221 . As we show in Appendix 
B, using the factorization approximation and an approx¬ 
imation for the steady state vertex function, we arrive at 
the following expression for the memory function, 


,,iriY + n P u ll(9) f dqidq 2 

M ( q;t ) = —— J 3 £(q —qi —q 2 ) (8) 

x (q- [qiC(gi) + q 2 C(q 2 )}) 2 F( qi ;t)F(q 2 ;t). 

Equation 0 has a structure similar to the memory func¬ 
tion of the mode-coupling theory, but it involves a new 
function C(q) (which replaces the direct correlation func¬ 
tion c(q) in the mode-coupling M 1TI (q;t)), 


pC{q) = 1 - 


^||(g) 

wy (00)5(9)’ 


(9) 


where W||(00) = (3N£q) 1 (^2 i (f) + F*) 2 ^. Equations 

00 are closed and can be solved if static steady state 
functions S(q) and wy (9) are available. To test the theory 
quantitatively, we used the static information obtained 
directly from simulations. Using the KA system would 
require formulating and solving our theory for a binary 
mixture. Instead, we performed additional simulations 
of a one-component Lennard-Jones (LJ) system of self- 
propelled particles [35]. We measured S(q) and wn ( 9) and 
used them to solve Eqs. 00 numerically to compare the 
dynamics predicted by the theory with the numerical re¬ 
sults. We focus on the observed non-monotonic evolution 
of the dynamics with the persistence time because it rep¬ 
resents a demanding test of the theory. 

In Fig0a-c) we show the dependence of g(r), S(q) and 
W|| {q)/uj\\ (00) of the self-propelled LJ system on the per¬ 
sistence time T p at constant T e ff = 0.9 (this is the lowest 
T e ff at which we were able to simulate our one-component 
system without observing spontaneous ordering). Again, 
the structure becomes progressively more pronounced as 
T p increases. Concurrently, correlations of particle veloc¬ 
ities also develop, as revealed by W||(9)/w||(oo). Finally, 
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FIG. 2: Structure (a,b) and dynamics (c) obtained from 
simulations, and predicted by the theory (d) for various t p 
at T e g = 0.9 in the one-component LJ system, (a) Pair dis¬ 
tribution function g(r ); the curves are shifted vertically for 
clarity; r p = 0, 2 x 10 —3 , and 2 x 10 -2 (bottom to top). In 
(b) thick lines represent wy (q )/wii (oo) and thin lines represent 
S(k) at t p = 2 x 10~ 4 (solid) and 2 x 10 -2 (dashed). The in¬ 
set in (b) shows the persistence time dependence of un (oo)r p . 
In (c) we show r a (circles, left axis) and the peak value of 
g(r) (triangles, right axis), obtained from simulations. In (d) 
we show r a (circles) and the inverse self-diffusion coefficient, 
0.1/.D (squares) predicted by the theory. 


the quantity wm (oo)r p , which is a measure of local mobil¬ 
ity in the interacting self-propelled system, decreases. 

Because it incorporates these different trends, our the¬ 
ory is able to account for the non-monotonic dependence 
of the dynamics on the persistence time. In particular, 
with increasing r p both S(q) and ui\\(q) grow for q around 
27r/cr and as a result C(q) gets smaller than c(q), which 
is the likely source of the speed-up in the dynamics. At 
larger r p decreasing W||(oo)t p and increasingly more pro¬ 
nounced local structure seem to prevail upon the increase 
in velocity correlations resulting in the slowing down in 
the dynamics. In panels (c,d), we show the dynamics pre¬ 
dicted by the theory and obtained from simulations, re¬ 
spectively. Clearly, our theory qualitatively predicts the 
non-monotonic dependence of the dynamics on the per¬ 
sistence time, suggesting that including non-equilibrium 
velocity correlations ui y ( q ) in the theory is of major im¬ 
portance. Mode-coupling theories overestimate the slow¬ 
ing down of the dynamics 122 ] , and this explains why our 
theory predicts a more pronounced non-monotonic effect 
than in the simulations. Describing quantitatively glassy 
dynamics at thermal equilibrium is a notoriously difficult 
and open challenge m- Therefore, quantitative agree¬ 
ment should not be expected in the far from equilibrium 
context of self-propelled particles. 

Various kinds of generalizations of mode-coupling the¬ 
ory for driven glassy fluids have been proposed nmm 


fIDl - fliT] . In particular, Ref. 03] developed a theory for ac¬ 
tive Brownian particles, where self-propulsion is added 
to a thermal Brownian system, but this work differs 
from our approach on important aspects. First, corre¬ 
lations between positions and self-propulsions were ne¬ 
glected in [23]. Technically, this amounts to replacing 
the local steady-state average in Eq. © by the aver¬ 
age over the distribution of self-propulsions. As a re¬ 
sult, ((Mj) lss - <f;) lss <fj) ta ) gets replaced by DfT p SijI , 
and the steady state correlation function W||(q), which 
we have shown to play a central role, does not appear. 
Second, in our derivation we use projection operators de¬ 
fined in terms of the steady state distribution, whereas 
Ref. [23] uses the equilibrium distribution. As a result, 
the memory function derived in Ref. [251 is the same as in 
the equilibrium mode-coupling theory while ours is dif¬ 
ferent. Because we consider an intrinsically nonequilib¬ 
rium system, there is no equilibrium distribution that we 
could use. Physically, the steady state distribution seems 
more natural since we are describing fluctuations in the 
steady state. An obvious disadvantage of our choice is 
that we have to obtain the steady state correlation func¬ 
tions, S(q) and ui\\(q), from simulations. 


V. SUMMARY 

Using computer simulations and tools from liquid state 
theory we developed and analyzed an athermal system 
of interacting self-propelled particles. We showed that 
the microscopic structure and long-time dynamics evolve 
nontrivially with the degree of nonequilibrium, which 
challenges equilibrium theories for dense fluids. We pre¬ 
sented a theory for the collective dynamics of an active 
many-body system that can qualitatively capture these 
phenomena. In particular, the speed up of the dynamics 
of the active fluid was linked to emerging steady state 
correlations of the velocities, an object with no relevant 
equilibrium counterpart. In future work, we will numeri¬ 
cally characterize the approach to the glass transition in 
more detail for such nonequilibrium conditions. On the 
theory side, we will analyze the nature of the ‘nonequilib¬ 
rium glass transition ’ m predicted by our theory. We 
shall also compare the dynamics predicted by the the¬ 
ory with that obtained from simulations at larger degrees 
of supercooling, and the relation between the ergodicity 
breaking temperature predicted by our theory with the 
transition temperature obtained from simulations. To 
perform detailed quantitative comparisons, we will need 
to generalize the present results to binary mixtures. More 
generally, our work paves the way for developing a gen¬ 
eral, microscopic understanding of the glassy dynamics 
of active materials when different interparticle interac¬ 
tions, self-propulsion mechanisms and possibly thermal 
noise are at play. 
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Appendix A: Derivation of effective evolution 
operator n eff (.z) 


where Pjv(ri,..., rjv; t) is the N-particle distribution of 
positions, 

Pjv( ri,...,rjv;t) = / dfi...dfjvPjv(iT, fi,.... rjv, f N -,t), 

(A5) 

and jj(ri,..., rjv; t) is the current density of particle i. 
j»( r i, •••, rjv; t) = (A6) 

to 1 J d fi ”.d f Ar (Fj + f i )Pjv(r 1 ,f 1 ,...,r JV ,fj V ;i)- 

Our assumption of the absence of systematic currents in 
the steady state is implemented as follows, 


In Appendix A we present an outline of a derivation of 
the effective IV-particle evolution operator (^r), Eq. 
([6]). There are three main assumptions used in this 
derivation. First, we assume that systematic currents 
vanish in the steady state of our system. Second, we as¬ 
sume a separation of timescales for the structural relax¬ 
ation and the relaxation of the self-propulsions. Third, 
we approximate the dynamics in the space orthogonal to 
the local equilibrium distribution of self-propulsions by 
the free relaxation of the self-propulsions. 

Equations of motion ([l]|2j) are equivalent to the follow¬ 
ing evolution equation for the A-particle joint distribu¬ 
tion of positions and self-propulsions, 


d t P N ( ri,fi, ...,r N ,f N ;t) = fiPjv(ri, fi, ...,r N ,f N ;t) 

(Al) 

where f l is the A-particle evolution operator given by 

Eq. 0. 

We assume that the evolution equation (Al) has a 
steady state solution Pjy (iq, f x ,..., r^v, fjv), and thus 

HP^ s (r 1 ,f 1 ,...,r JV ,fjv) = 0. (A2) 


In the main text and in the following we use brack¬ 
ets (...) to denote averaging over the joint steady 
state distribution of positions and self-propulsions, 
Pjv ( r l) fi, rjv, fjv)- 

From the joint steady state distribution we can obtain 
a steady state distribution of positions, Pjy(iq,..., rjy), 


P/v( r i, •••, r jv) = I dii...di N P^(r 1 ,{ 1 , 

(A3) 

In the main text and in the following we use brackets 
(... ) r to denote averaging over a steady state distribution 
of positions, Pj^( r 1 ,...,r N ). 


We assume that in the steady state there are no sys¬ 
tematic currents. To make this statement more precise 
we first define the current density by integrating Eq. (Al) 
over self-propulsions to get the following continuity equa¬ 
tion, 


d t Pjv(ri,...,rjv;i) = V; • ji(iT,...,r^; t), (A4) 


AT 1 / ,/f, c/l\ IF - f P??fr 


The above equality implies that the local steady state 
average of the self-propulsion is equal to the negative 
of the force, (f*) lss = — Fj, where the local steady state 
average is defined as 

(••■)i ss = (A8) 

We assume that the self-propulsions relax faster than 
the positions of the particles. This assumption is appli¬ 
cable for strongly interacting systems where structural 
relaxation is slowing down, whereas the evolution of the 
self-propulsions stays, by definition, independent of inter- 
molecular interactions. The separation of the timescales 
for the structural and self-propulsions relaxations allows 
us to derive an approximate equation of motion for the 
distribution of particle positions. 

We define the projection operator on a local 
equilibrium-like distribution ( i.e . on a distribution in 
which self-propulsions have a steady-state distribution 
for a given sample of positions), 


PissF/v(ri,fi,...,rjv,fjv;i) 
= Pf (ri,fi,...,rjv,fjy) 
Pf (ri,..., r/v) 

= P/v ( r i, fi, ■ ■ ■, r jy, fjy) 
Pjv (d, ...,r N ) 


J dfi...dfjvPjv(ri, fi, r N , f N ; t) 

Pjv(ri, •••, rjv; t). (A9) 


Next, we define the orthogonal projection, 
Qiss = X — Piss, and write down equations 
of motion for PissPv (r-| , f-] ,..., r ,v, f/v ; t) and 
Qis S P/v(ri,fi, ...,rjv,f N ;t), 


9tPi ss P/v(ri, fi,..., rjv, fjv; t) = 

PissHPissPvjri, fi, r N , fjv; t) 
+Pi ss HQi ss Pjv(ri,fi,..., rjv, f)v; t), (A10) 


9tQi ss Pjv(ri,fi, ...,r N ,f N ;t) = 
QissHPis S Piv(ri, fi,..., rjv, fjv; t) 
XQlss^QlssPN (^T; fl? ■ • ■ ■ I*jv, f N] 1) ■ 


(All) 









6 


Since our final goal is to calculate the intermediate 
scattering function given by Eq. ([3]), we can assume 
that Qi ss -P/v(r 1 , fi,rjv, fjv; t = 0) = 0. Then we can 


solve Eqs. ( A10|All ) for the Laplace transform, £T, of 
9*Piss-Pw(rii fij r N , t), which is given by 


CT [d t Vi ss P N (ri,fi, ...,r N , f N ;t)] (z) = 


Plss^Plss + Piss ^ Sis 


z Qiss^Qis 


-Qlss^Pls; 


P ls sPjv(ri,fi, ■ rjv, fjv! z). 

(A12) 


Using the assumption that systematic currents vanish in the steady state, Eq. (A7), one can show that the first term 

(A13) 


inside the brackets on right-hand-side of Eq. (A12) vanishes. Furthermore, one can show that 

Pjv(ri, ...,r N ]z) 


SissfiPissPv(^) = -4 _1 E( f *- < f i)u») iJ v(r 1 ,f 1 ,-,r JV ,fj V ) 


_ 1 P|f(ri,..., rjv) . 


and 


PlssUQlss- = - P "pI/^Z7n) N) ^ ^ Vi ' / dfl - dfN - 


(A14) 


We note that Q\ ss FlQi S s describes evolution in the 
space orthogonal to the local steady state space. We 
assume that this evolution is entirely due to the 
free relaxation of the self-propulsions. Specifically, 
we assume that Qi ss UQi ss can be approximated by 
( t~ x U + We note that this approxi¬ 


mation physically means that the relaxation rate of the 
self-propulsions is not renormalized by the interparti¬ 
cle interactions. Combining the last approximation with 
Eqs. ( A13]|A14 | we get the following approximate equal¬ 
ity 


PissUQi ss (z — Qi S sUQi ss ) 1 Qi ss nP lss P/v(iT,fi, ...,r N ,f N ;z) « 

P“(ri,fi)- ! riv,fjv) t _2 v fir ,c / f \ x 

^(f^(fi) lss )PA s (ri,fi,...,riv,fiv)- 

l 

-1 

and integrate by parts. Finally, we integrate both sides of the 
resulting equation over self-propulsions and get the following expression for the Laplace transform of dt-Pjv( r i,r a/-; t) 

CT[d t P N {v u ...,v N -t)} (z) =4“ 2 E V * • (z + Tp 1 )' 1 (( f i f t>i ss - (fi)iss ft >1 J • [—Ff + V,] P N (r 1 ,...,r N -,z){A16) 


Now, we expand 


N 

z - £i=i 


JL 

dfi 


D 


f dU 


P N (ri,...,r N ;z) 
1 P s J{r i,...,rjv) 


. (A15) 



The right-hand-side of Eq. (A16) defines the effective 
evolution operator fI off (;z) given by Eq. 


Appendix B: Derivation of an approximate evolution 
equation for F(q\t ) 


In Appendix B we present an outline of a derivation of 
an approximate evolution equation for the intermediate 
scattering function F(q;t ), Eqs. (7|9). The framework 
of the derivation is based on that of the derivation of 
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the mode-coupling theory: first, the Laplace transform 
of the time derivative of F(q;t) is expressed in terms of 
the frequency matrix and the reducible maun memory 
matrix. Next, the reducible memory matrix is expressed 
in terms of the irreducible one. Finally, an approximate 
expression for the irreducible memory matrix in terms of 
the intermediate scattering functions is derived. 

The Laplace transform of the intermediate scattering 
function defined in Eq. (|3j) reads 

F(q-z) = N^'(n{q)(z-Cl)~ 1 n(-q)^ 

= N~ 1 (n(q) (z - n(-q)) r (Bl) 

where n(q) = JT e _?q ' r ' is the Fourier transform of the 
microscopic density and in the second equality we used 
the effective evolution operator derived in the previous 
section of this supplementary material. 

To derive an approximate evolution equation for in¬ 


termediate scattering function F(q;t) we first define a 
projection operator on the microscopic density 

V n = ... n(-q)) r (n(q)n(-q)) I T 1 (n(q)... . (B2) 

We should emphasize that projection operator V n is de¬ 
fined in terms of the steady state distribution, unlike in 
the approach of Farage and Brader [23]. Next, we use 
the identity 


z-VL eS {z) z-£l eS {z)Q n 

_l _ 1 _ Q, eS ( z )V _ 1 _ 

z-n° s ( z )Q n z-n eS (z)’ 

where Q n is the projection on the space orthogonal to 
that spanned by the microscopic density, Q n = I — V n , 
to rewrite the Laplace transform of the time derivative 
of NF(q;t) as follows 


CT[d t NF(q;t)}(z) = (n(q)n ett (z) 


z — 


(-q) = (n(q)Q, ett (z)V., 


'z- 


i(-q) 


n(q)n cS (z)Q n ^ _ ^ cff ^ n(-q)^) = ( rl (q)^ eff (-)^(^q)) r (n(q)n(-q)) r 1 (n(q) ^ n(-q) 

1 


n(q)fl cS (z)Q 


Q n tt cS (z)n{- q) ) (n(q)n(—q)) r 1 ( n(q) 


-1 


1 z - Q n n^(z)Q, 

We identify the analogue of the frequency matrix, z). 

q • (E,, m )i ss - (fi)iss J e-^-^ \ ■ 


z — 


n(-q) 


(B4) 


<n(q)n cff (2:)n(-q)) r = -- 


Co (~ + t p *) 


= ~q 2 N = -q 2 NH(q-,z)( B5) 

[Z + Tp ) 


where W||(g) is the function quantifying correlations of Note that ((■•■)i ss ) r = (■••)• Furthermore, we identify 
the velocities of individual particles that was introduced the analogue of the reducible uni HU memory matrix, 


below Eq. (J7|, 
wy(g) = 


M{q;z), 


^ • (^2 ((f,f,) lss - (f,) lss (f,) J e-^-* 




= q • (J2 & + F ») ( f i + F ^) e_iq ' (r<_ri) y • q- (B6) 


n(q)n cS (z)Qn z _ Q?i ^eff (3 )Q n Q»n Gff (^)n(-q)^ = (Co (* + r p E) q ■ (y, e iqr ‘ (( f * f J>i ss “ ( f *>iss &K) 


x [-V,+Ff] Q„ — 


^ eS{z)Q Qn E • « f ^>i ss - < f *Es (WiJ ) • q = q 2 NM(q ; z) 

1 ^ n l,m / r 


(B7) 





















Following the procedure used previously to derive the 
mode-coupling theory for Brownian systems we introduce 
irreducible memory matrix A4 lrl (q;z), which is given by 
the expression analogous to Eq. ( |B7[ ) but with the pro¬ 
jected evolution operator Q n £l eS (z)Q n replaced by irre¬ 
ducible evolution operator 9 lrr (;r) j40M42] . The relation 
between M{q\z) and M m (q\z) reads 


M{q; z) = A4 irT (q; z) - M m {q ; z)n~\q] z)M{q ; z). 

(B8) 

Combining Eqs. (B4 B5) and ( |B7]|B8 ) we can write 
the Laplace transform of the time derivative of the inter¬ 
mediate scattering function in the following way, 


CT[d t F{q-t)\{z) = (B9) 

- q 2 n(q ; z) (1 + M irr (q; z)/H(q ; z))^ F(q ; z)/S(q) 

where S(q) is the steady st ate s tructure factor, S(q) = 
(n(— q)n(q)) r . In turn, Eq. (B9) can be re-written as 


(z + Tp 1 + (z + T p x ) 2 M iri (q;z)/uj ||(g)) (BIO) 

X ( zF(q;z ) - F(q;t = 0)) = - (wy (q)q 2 /S(q)) F{q;z). 


Eq. ( |B10| transformed back to the time domain repro¬ 
duces Eq. |7]). In particular, the inverse Laplace trans¬ 
form of (z + Tp 1 ) M 1IT (q;z)/uj\\(q) is equal to the irre¬ 
ducible memory function M m (q\t) that enters into Eq. 

0 - 


To proceed we derive an approximate expression for 
the irreducible memory function M' rv (q\ t) in terms of the 
intermediate scattering functions. To this end we first fol¬ 
low the steps of the derivation of the mode-coupling the¬ 
ory m and replace projection operators Q n in M lrr (g; t) 
by projections on density pairs, and next factorize both 
steady state and time-dependent four-point correlations. 
As a result we get the following approximate expression 
for M' lr (q- : t), 


M^(q-t) 




iVy E q • (e ■ e_iq ‘ ri (&Viss - <f i)iss &•>«„) • [-Vi + Ff] Q„n 2 (- qi , -q 2 )\ 

1 Id 1 qi,...,q s \ ij / 


5qiq35q2q4 + S « lClJ ™ 3 -N 2 F(q 3 -,t)F(q,-,t) (6 a . a .S, 


qaqs^qjqe + ^qe^nqs 


. ^q5q7^' 


N^S( qi )S(q 2 ) 

'N^qjslq^ ( n2(q7 ’ q8)S " ^ Vfc ‘ mi)lss ~ {ik)lss 6iqr 


q 


(Bll) 


Next, we approximate the vertex functions. The justi- discussed elsewhere. Here we present the approximate 
fication for the form of this last approximation will be expression for the left vertex, 


C 2 q • (E' e_<q ‘ ri - &)us ft), J ' [~Vi + Ff ] Q„n 2 (-q 1; -q 2 )\ 

\ / r 


-iqqiN \^Ji\(q) ^ i\{q 1 )S{q 2 ) - Wy (q)S(q 1 )S{q 2 ) j ^q,q 1+ q 2 

~iq q 2 N (wy(g)—^-xWy (<? 2 )S(<?i) - wy (q)S(q 1 )S(q 2 )) 6 q , qi+q2 , 
V wy(oo) J 


(B12) 


■ L g—>00 Wy ( q) 


where wy ( 00 ) = lim, 

Finally, we use Eq. (19]) and an analogous approxima- 


= tion for the right vertex in Eq. (Bll) and, after taking 


the thermodynamic limit, we obtain Eqs. (|8ll9|). 
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